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I. INTRODUCTION 



High-precision measurements of the ground-state hyperfine sphtting (HFS) were per- 
brmed for various heavy H-hke ions, including ^°^Bi, ^^^Ho, ^^^Re, ^°''Pb, ^osrjj g^^^ 205rp| 
5|. Progress in experiments motivated intensive theoretical calculations of the hyper- 
fine splitting in highly charged heavy ions aiming to test quantum electrodynamics 
(QED) in strong electromagnetic fields. It was found that in heavy ions the QED effects 
are obscured by the uncertainty of the nuclear magnetization distribution correction (Bohr- 
Weisskopf effect). However, simultaneous study of H- and Li-like ions of the same isotope 
can help to overcome this problem, since the uncertainty of the Bohr- Weisskopf effect is 
significantly reduced in the specific difference of the corresponding HFS values 13]. High- 
precision measurements of the hyperfine splitting of H- and Li-like bismuth are feasible at 
the experimental storage ring (ESR) and the HITRAP facility in GSI {20! . 21 1. Recently, 
after 13 years of attempts, the HFS of the ground state Li-like Bi has been directly observed 



in GSI 



22| . These measurements together with accurate theoretical calculations will provide 



the possibility for the stringent tests of QED in strong fields. 

The recent improvements of the theoretical accuracy for the specific difference of the 
HFS values are related to the calculations of the screened QED corrections 23, 241 and the 

n 

two-photon exchange corrections 25| to the HFS of the Li-like ions. Now the uncertainty 



of the specific HFS difference is main 



screened QED corrections. In Refs. 



16 



y determined by the Wichmann-Kroll part of the 



26143011 the screened QED corrections were evaluated 



by introducing an effective local screening potential in the zeroth-order (Dirac) equation. 
However, the screening potential approximation does not provide reliable estimation of the 
uncertainty. Recently, the two-electron self-energy diagrams and a dominant part of the 
two-electron vacuum-polarization diagrams, which represent the leading contribution to this 



effect, have been evaluated within the systematic QED approach 



23 



2J] . The present paper 



ectron vacuum-polarization diagrams. 



24| . In particular, the electric-loop and 



is devoted to the rigorous evaluation of the two-el 
which have been treated approximately in Refs. 23, 
magnetic-loop diagrams are evaluated to all orders in aZ, including the Wichmann-Kroll 
contributions. The Wichmann-Kroll terms of the remaining internal-loop contributions are 
estimated. The numerical results are presented for the hyperfine structure of Li-like bismuth 

209gj80+ 
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The relativistic units {h = 1, c = 1, m = 1) and the Heaviside charge unit 
[a = e^/(47r), e < 0] are used throughout the paper. 



II. BASIC FORMULAS 

The interaction of atomic electrons with the nuclear magnetic moment is described by 
the Fermi-Breit operator, 



i^. = ^M-T, (1) 



where ^t is the operator of the nuclear magnetic moment. The electronic operator T is given 
by 

T = E -^^nr^) , (2) 



where the summation runs over the atomic electrons, ex. is the Dirac- matrix vector, rij = rj/r. 



and F{r) is the nuclear magnetization distribution factor (see, e.g., Refs. [Ul, |30|). Here we 
employ the homogeneous sphere model, 

1 /2 

where Rq is the radius of the sphere, related to the root-mean-square charge radius (r^) of 
the nucleus as Rq = A/5/3(r^)^^^. The ground-state hyperfine splitting of a highly charged 
Li-like ion in the non-recoil limit can be written as 



(„) _ ajaZY gj 



A{aZ){l - -e) + ^B{aZ) + ^C{aZ) 
+ ^D{Z, aZ) + xqed + a;sQED 



(4) 



Here gi = ^/{^nI) is the g factor of the nucleus with magnetic moment /x and spin /, 
/iAT is the nuclear magneton, and nip denotes the proton mass. A{aZ) is the one-electron 
relativistic factor, 6 and e are the corrections due to the finite distribution of the charge and 
the magnetic moment over the nucleus, respectively, which can be found either analytically 



m 



\3l\ or numerically. The interelectronic-interaction correction of first order in 1/Z is 
represented by the function B{aZ). The function C{aZ) incorporates the interelectronic- 
interaction corrections in the second order in 1/Z, D{Z,aZ) corresponds to the third- and 
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FIG. 1: Diagrams, contributing to the screened VP correction to the hyperfine sphtting. The wavy 
Une indicates the photon propagator and the double Une indicates the electron propagator in the 
Coulomb field. The dashed line terminated by the triangle denotes the interaction with the nuclear 
magnetic field. 

higher-order corrections in 1/Z. xqed and xsqed correspond to the one-electron and many- 
electron (screened) QED corrections, respectively. 

To first order in a and 1/Z, the screened QED correction xsqed to the hyperfine splitting 
is given by the sum of the self-energy (SE) and vacuum-polarization (VP) parts 



SE VP 
3:^SQED = 3^SQED + ^SQED" 



(5) 



In the present paper the vacuum-polarization part Xsqe^i is considered. The corresponding 
diagrams are depicted in Fig. [TJ The total contribution of these diagrams is conveniently 
divided into reducible and irreducible parts. The irreducible parts of each diagram A-F 
are denoted by the same letter: a;gQgQ'^''. The reducible contributions are to be considered 
together with the non-diagram terms (see Ref. 2^ for details). They are divided into three 
parts, G, H, and I, according to the type of the vacuum-polarization loop: G is associated 
with the electric-loop diagrams A, B, and E; H corresponds to the magnetic-loop diagram 
C; and I — to the internal-loop diagram F. 
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The total correction due to the two-electron vacuum-polarization is given by the sum 



™VP _ ^VP(A) VP(B) VP(C) VP(D) VP(E) VP(F) 

•^SQED — "''SQED '''SQED "''SQED "''SQED "''SQED "''SQED 

, VP(G) VP(H) VP(I) 

^"''SQED ' "''SQED ' "''SQED • \^) 



A. Electric-loop diagrams 



The terms A, B, E, and G in Eq. (jH]) correspond to the so-called electric-loop (el) diagrams 
A, B and E. The contributions of these terms are given by the following expressions 



X 



VP(A) 
SQED 



b PQ 



ni,n2 ^ 



{n^\To\n2){n2Pb\I{A)\QaQb) 



„i , .{niPb\I{A)\n2Qb){n2\To\Qa) 



(7) 



aQb) 



b PQ 



nx,n2 



£Pa — ^rii 



^Pa ^nn 



(8) 



b PQ 



Y^\Pa\UUn,) 



n\,n2 



el ,^ , {Pb\ To\n2){n^n2\I{A)\QaQb) 

£Pa — £ni){£Pb — £712) 

{n,Pb\I{A)\Qan2){n2\To\Qb) 



+ > {Pa\Uyp\ni) — 



){eQb - ^nj 



, (9) 
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VP{G) 
"''SQED 



b PQ ni 

{ (Pa|f/^'pl^i) (^1 l^ol^'a) {PaPh\I{/\)\QaQh) 
+ {Pa\U^^\Pa){Pa\TQ\ni){niPb\I{A)\QaQb) 
+ {Pa\TQ\Pa){Pa\U^^\ni){niPh\I{^)\QaQh)] 

b PQ I ni ^P'' 

X [{Pa\U^j.\ni){ni\TQ\Pa){PaPb\I\^)\QaQh) 

+ {Pa\U'^p\n^){nrPh\I'{I\)\QaQh) {{Qb\To\Qb) - {Pb\To\P 
{PaPb\riA)\niQb){ni\To\Qa) 



+ {Pa\U^^p\Pa) 



{PaPb\riA)\Qani){ni\To\Qb) 



EQb - En, 
b PQ 

X {{Qb\To\Qb) - {Pb\To\Pb)) . (10) 

Here and below (Eqs. f lTB]) . f lT7|) . f l5U]) . (EI])) -P and Q are the permutation operators, 
interchanging a and 6, A = eg;, — ^qq. The summation over 6 runs over two core electron 
states with different projections of the angular momentum. The prime at the summation sign 
Y^' indicates that the energy of the intermediate state differs from the energy of the initial 
state, so that the corresponding denominator is non-zero. The interelectronic-interaction 



operator /(A) and its derivatives are defined as in Ref. 32] • The factor Ga is defined by 

the quantum numbers of the valence state, 

_ n^(2/ + l)j(j + l) 3 

2{aZfm, {aZfrrij' ^ ' 

where irij is the projection of the angular momentum j. 

Equations fl7|)-f lT0|) involve the matrix elements of the standard electric-field-induced 

vacuum-polarization potential U^-p- The unrenormalized expression for U^-p is given by 

UU^) = 7^ [ ^'r'^- r duTr[Giuj,r',r')] , (12) 
where G{u, r, r') is the Dirac-Coulomb Green function 
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FIG. 2: The decomposition of the C/yp into the Uehhng (Ue) and the Wichmann-Kroll (WK) 
terms. 



The decomposition of the vacuum-polarization loop into the Uehling (Ue) and the 
Wichmann-Kroll (WK) terms is depicted in Fig. [2l In this expansion only the lowest order 
term, the Uehling term, is divergent. The charge renormalization yields a finite well known 
renormalized expression 



el-Ue/ 
VP I 



-aZ — / dr' Anr' p(r' 
1 



X— [exp(— 2|r — r'\t) — exp(— 2(r + r')t)] 



(14) 



where the density of the nuclear charge distribution p(r) is normalized to 1. The Wichmann- 
Kroll part (the term in brackets in the decomposition in Fig. |2]) is calculated to all orders 
in aZ as the difference between the unrenormalized total and Uehling contributions. The 
resulting expression is free from divergencies, completely isolated in the Uehling term 33 1. 
It is known, that no spurious terms contribute, if the calculation is based on the the partial- 
wave expansion of the electron Green function and the summation is terminated after a finite 



number of terms 



34 



36l ]. Thus the WK contribution after Wick- rotation of the contour of 



w-integration in the complex plane can be written as 



37| 



[/^M = ^Vkl / du I dr'r" 



1 II Ii2 

dr r - 



maxi^r, r 



)]} 



(15) 



,k=l 



Here Gjf and F**"' are the radial components of the partial-wave contributions to the bound 
and free electron Green's functions, respectively, and V{r) is the electric potential of the 
extended nucleus. 
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B. Magnetic- loop diagrams 



by 



he magnetic-loop (ml) diagram C corresponds to the terms C and H in Eq. (I6j), given 



23 



24| 



VP(C) 
"''SQED 



b PQ ni ^ 



b PQ 



Equations f|T6|) . f|T71) involve the matrix elements of the magnetic- field-induced vacuum- 
polarization potential Uyp. Its unrenormalized expression reads 

f^vp(r) = ^ r duj [ dh' [ d3r"_^Tr[aG(u;, r', v"W)G{u, r", r')] • (18) 
2vr2 J J |r - r'l 

The scalar product is implicit in Eq. f|T8|) . Similar to the electric- loop potential, the decom- 
position of the magnetic-loop potential into the first order term and the remainder (Fig. 
E]) leads to the isolation of the divergency in the leading Uehling term. It is given by the 
Eq. ( ITSl) with the bound-electron Green function replaced by the free-electron one. For the 
sphere model of the nuclear magnetization distribution {F{r) given by Eq. ([3])) the analytical 



expression for the renormalized magnetic-loop Uehling term reads 
rml-Ue^„^ _ a [n X q:]o 3_ 

'^0 



3Q| 



U^p{r) = ^777753 [4ri?o[/3i(i?o + r) + /3i(|i?o - r\ 



+2{Ro + r)l32{Ro + r)- 2\Ro - r|/32(|i?o - r|) 
+f3s{Ro + r)-/3s{\Ro-r\)] , (19) 



where the function /3„ is defined by 



/3n(r) = 1 Tdt f 1 + 7^ ) exp(-2tr) . (20) 



3 ii t"+2 V 2t2^ 

Similar to the case of electric-loop term, the magnetic-loop Wichmann-Kroll contribution is 
calculated by summing up the partial-wave differences between the unrenormalized total and 
Uehling contributions. The magnetic-loop diagram contributes also to the nuclear magnetic 
moment. The corresponding Uehling term is equal to zero, however the Wichmann-Kroll 
term is not. Therefore, in the calculations of the WK contribution one should account for 
the related contribution to the nuclear magnetic moment in the zeroth-order HFS value (see 
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Refs. 12|, |l8|, |38|). This implies the replacement of the nuclear magnetic moment fi by the 
'bare' value 

/i /ibare = /i - A/i . (21) 

The nuclear magnetic moment correction A/i due to the magnetic-loop WK part can be 
expressed as Afi = efi with the dimensionless parameter e given by 
1 a 



2m 2 



/p poo 
dh / d^r' / du [Tr {[r x a]oG(w, r, r')ro(r')G(w, r', r)} 

-Tr{[r X a]oF(a;,r,r')To(r')F(a;,r',r)}] . (22) 



Finally, the corrected magnetic-loop WK contribution is obtained by subtraction of cTq from 
the WK part of Eq. flTSl) . The corresponding expression reads 



f/^™(r) 



a 

X 



— oo 

a 



doo I dh' / d-'t 



-[Tr{«G(c^,r',r")To(r")G(a;,r",r')} 
r — r'l 

-Tr{aF(a;,r',r")ro(r")F(u;,r",r')}] 

-iTo(r) [Tr{[r' x a]oG(a;, r', r")To(r")G(a;, r", r')} 

-Tr {[r' X a]oF(u;, r', r")To(r")F(a;, r", r')}] ) • (23) 

The angular-momentum conservation allows one to distinguish in this expression the follow- 
ing angular dependence 

f/^™(r) = [n X a]o <™(r) . (24) 

The radial part MYp^^(r) can be presented in the form 

<™(r) = dr'r" - ^ F{r)^ p^,{r') , (25) 

where the magnetic-loop WK charge density PnixiT') after rotating the contour of uj- 
integration can be written as 



poo poo 

^ dr"F{r") duTB'iK^,K,)S,,,,iu,r',r"). 

6^ Jo Jo rf„ 



(26) 



Here the Green-function trace Skih2 is given by 

5.,«,(u;,r',r") = Re {Gi\{zu,r' ,r'') Gll{zu,r' ,r'') + Gll{zu,r\r'') Gl{zu;,r' ,r'') 

+Gl\itu, r', r") Gllitco, r', r") + G^^^{iuj, r', r") G'^^iiu:, r', r") 



I "\ 77122/ • / //\ 77112/- / r/\ 77i21/- / r/\ 



F-(za;,r',r")F-(^a;,r',r")K2r) 




o 



< + 




FIG. 3: The decomposition of the C/yp into the UehUng and the Wichmann-Kroh terms, 
and the angular factor B{k,i,K2) is defined by the following expression 

B{^„^2) = —^ v/2(2ji + l)(2j2 + l)(-l)M _^ 1 . (28) 

The selection rule, following from this expression, forces k,2 to be equal to either ki or — ki±1. 
Therefore the double ^-summation can be reduced to a single one, including diagonal and 
off-diagonal K-terms, 



where k = —k, — sienf/t). 



- 1/4 



+ 1/2 



(29) 



C. Internal-loop diagrams 



The internal- loop diagram F corresponds to the terms F and I in Eq. ([6]), given by 231.124 1 



= 2G.$:$:(-l)--«$:'i^^^(n,P5|/vp(A)|gaQ6) , (30) 



b PQ 



Sd = ^G„5^5^(-l)^+«(PaP6|J(.p(A)|QaQ6) {{Qb\To\Qb) - {Pb\To\Pb)) . (31) 

b PQ 

Equations (I30p and (EI]) involve the interelectronic-interaction operator, modified by the 
vacuum-polarization loop, 

> 2 /"OO 



/vp(e,ri,r2) = ^ / du I dh[ / d% 



ai^exp(z 


e\ 




- r'J) a2„ex.p{i 




|r2 


-r'^l) 












r2 




r'2 





xTr [a^Giu - e/2, r[, r^)a''G(a; + e/2, r^, r'^)] 



(32) 



where e is the energy of the transmitted photon, = {I, ex.) is the four- vector of the Dirac 
matrices, and the summation over fi and u is implicit. The matrices and act on 
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FIG. 4: The decomposition of the internal electric- loop /yp into the Uehling and the Wichmann- 
Kroll terms. 

the spinor variables corresponding to ri and r2, respectively. This operator is also divided 
into the leading divergent Uehling part and the remaining finite Wichmann-Kroll part (see 
Fig. H]). The renormalized expression for the Uehling term reads (see, e.g., js^l) 



dt 1 



1 

2t2 



exp -V4t2 - £2 In - ral . (33) 



37r |ri - rsl Ji V' ' J t"^ 
The corresponding contribution is taken into account rigorously. The direct part of the 
Wichmann-Kroll contribution with the hyperfine interaction vertex on the 2s-electron line 
has been calculated by introducing the screening potential of the (Is)^ closed shell electrons 
into the electric loop. We take the difference between the contributions with the Green 
function inside the loop calculated with and without the screening potential. The calculation 
of the direct part of the diagram D, which involves the internal loop modified by the hyperfine 
interaction vertex, has been performed in a similar way. It is evaluated by including the 
potential of the closed shell into the Green functions of the internal loop, that leads to a 
diagram of the magnetic-loop type. The evaluation of the remaining exchange diagrams is 
currently underway. 



III. RESULTS AND DISCUSSION 



The contribution of the two-electron vacuum-polarization diagrams to the hyperfine split- 
ting of Li-like ion is calculated in coordinate space according to the formulas presented above. 
The electric-loop terms (A, B, E, G) and the magnetic-loop terms (C and H) are taken into 
account completely, including the Uehling and the Wichmann-Kroll parts. The internal-loop 
diagram F and the corresponding reducible term I are calculated rigorously in the Uehling 
approximation. The related WK contributions are evaluated for the direct parts with the hy- 
perfine interaction vertex on the 2s-electron. The remaining terms are estimated employing 

11 



the assumption that the ratio of the WK part to the corresponding Uehhng part is the same. 
The uncertainty is ascribed as high as 200% of the value obtained in this way. The contribu- 
tion of the direct part of the diagram D is evaluated as described in the previous section. We 
ascribe a 200% uncertainty to the D contribution as a conservative estimation of the uncal- 
culated exchange term. The numerical evaluation of the one-electron wave functions for the 
initial (a) and intermediate (7^1,2) states is performed using the dual kinetic balance (DKB) 



approach 



401. The Fermi model for 



with the basis set constructed from the B-splines 
the nuclear charge distribution is employed in these calculations. The Uehling parts of the 
vacuum-polarization potentials are calculated according to the expressions ffH]) . ffH?]) . and 
f l33p . The Wichmann-Kroll parts involve the free- and bound-electron Green functions for 
the Dirac equation. The spherical shell model (possessing analytical solution for the bare 
nucleus case) and Fermi model of the nuclear charge distribution were employed in calcu- 
lations of the Green function. The nuclear magnetization distribution effect is taken into 
account within the homogeneous sphere model. The partial- wave expansion is terminated 
at = 5 in case of the electric-loop Wichmann-Kroll parts and at |k| = 10 in case of the 
magnetic-loop Wichmann-Kroll parts. The remainders of the k-summations are estimated 
using the least-square inverse-polynomial fitting. The numerical calculation procedure has 
been performed in different gauges, the gauge invariance should hold for the complete set of 
the two-electron vacuum-polarization diagrams. The total values obtained in the Feynman 
and Coulomb gauges for the photon propagator mediating the interelectronic interaction 
agree within the level of the numerical accuracy. The results for Li-like bismuth 209gj80+ 
both gauges are presented term by term in Table HI 



In Table Hn the individual terms are compared with those from Ref. 2J]. For the electric- 
loop terms the sums of the Uehling and Wichmann-Kroll contributions are presented. In 
Ref. ||2^] the Wichmann-Kroll parts were calculated by means of the approximate formulas 
for the electric WK-potential from Ref. j4l|. For the magnetic- loop terms (C and H) the 
Wichmann-Kroll part is figured out separately (WK-ml) in order to compare the result with 
the estimation given in Ref. 2S[- In that work it was evaluated utilizing the hydrogenic 2s 



1^ 



assuming that it enters with the same screening ratio as the Uehling 



value from Ref. 
term. 

The uncertainty of the present evaluation is determined by the uncalculated terms in sets 
D and F. It was estimated as high as 200% of the calculated WK terms in these sets. 
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TABLE I: Screened VP corrections to the HFS of Li-like bismuth 209gj80+ terms of xsqed for 
the Coulomb and Feynman gauges. Ue = Uehling, WK = Wichmann-Kroll. 







Feynman 




Coulomb 




Ue 


WK 


Ue 


WK 


A 


-0.000 


5074 


0.000 


0178 


-0.000 5085 


0.000 0179 


B 


-0.000 


2192 


0.000 


0056 


-0.000 2166 


0.000 0055 


C 


-0.000 


1692 


0.000 


0461 


-0.000 1670 


0.000 0455 


D 






0.000 


0021 




0.000 0021 


E 


-0.000 


0033 


0.000 


0002 


-0.000 0031 


0.000 0002 


F 


0.000 


0015 


-0.000 


0002 


0.000 0015 


-0.000 0002 


G 


0.000 


2896 


-0.000 


0123 


0.000 2879 


-0.000 0124 


H 


0.000 


0023 


-0.000 


0006 


0.000 0001 


-0.000 0000 


I 


0.000 


0000 


0.000 


0000 


0.000 0000 


0.000 0000 


Total (A-I) 


-0.000 6056 


0.000 0586 


-0.000 6056 


0.000 0586 


Total (Ue+WK) 




-0.000 5470 




-0.000 5470 



In Table ITTTl the specific difference of the ground-state hyperfine splitting in H-like bismuth 
209Bj82+ Lj_ijj^g bismuth 2093^80+^ ^/^ ^ ^^(2s) _ ^/\^e^is\ is considered. The 

parameter ^ is chosen to cancel the Bohr-Weisskopf correction, ^ = 0.16886 j^]- The 
rms radius was taken to be (r^)^/^ = 5.5211 fm 42|, the nuclear spin and parity I'^ = 
9/2—, and the magnetic moment fi/fiN = 4.1106(2) j43|. The most accurate values for 

251, where the 



the interelectronic-interaction contributions are taken from the recent paper 
contribution of the two-photon-exchange diagrams has been evaluated in the framework of 
QED. The contribution of the screened vacuum polarization calculated in this work equals 
—0.188(2) meV that agrees with the previous value —0.187(6) meV from Refs. 23|, |24]. The 
uncertainty of the present result, being 3 times smaller than the previous one, is determined 
by the conservative estimates for the contributions which have not been taken into account 
rigorously so far. The first error bar in the total value of the specific difference —61.320(4) (5) 
originates from the uncertainties of the screened VP contribution and the 1/Z^ and h.o. 
interelectronic-interaction term. The second uncertainty comes from the nuclear magnetic 
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TABLE II: Screened VP corrections to the HFS of Li-like bismuth 209-gj80+ terms of xsqed in 



2j]. TheWichmann- 



the Feynman gauge are compared to the corresponding corrections from Ref. 
Kroll part of the magnetic-loop diagram (WK-ml) is separated for comparison with the estimation 



of Ref. 



24|. 



This work Ref. [24] 



A 


-0.000 


4897 


-0.000 4881 


B 


-0.000 


2136 


-0.000 2128 


C 


-0.000 


1692 


-0.000 1691 


D 


0.000 


0021(42) 




E 


-0.000 


0031 


-0.000 0031 


F 


0.000 


0013(3) 


0.000 0015 


G 


0.000 


2773 


0.000 2766 


H 


0.000 


0023 


0.000 0023 


I 


0.000 


0000 


0.000 0000 


WK-ml 


0.000 


0454 


0.000 05(2) 


Total 


-0.000 547(4) 


-0.000 54(2) 



4J|, and other nuclear effects, which are not 



moment, the nuclear polarization corrections 
completely cancelled in the specific difference. 

In summary, calculations of the major part of the screened vacuum-polarization correction 
to the hyperfine splitting in Li-like bismuth have been performed. The Wichmann-KroU 
contributions to the two-electron electric-loop and magnetic-loop diagrams and to the direct 
parts of the internal-loop diagrams have been evaluated within the rigorous QED approach. 
As a result, the accuracy of the screened QED contribution to the hyperfine splitting of 
Li-like bismuth has been significantly improved. These results, combined with the recent 
rigorous calculations of the two-photon exchange contributions 251], provide a new value of 



the specific difference of the HFS values in H-like and Li- 
of magnitude more precise compared to the previous one 



ike bismuth which is by an order 



23 



2M- 
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TABLE III: Individual contributions to the specific difference A.' E of the hyperfine sphttings of 
Li-hke and H-like bismuth ^"^^Bi. Units are meV. 



Dirac value 844.829 876.638 -31.809 

Interelectronic interaction, ~ 1/Z — 29.995 —29.995 

Interelectronic interaction, ~ 1/7? 0.258 0.258 

Interelectronic interaction, ~ 1/Z'^ and h.o. — 0.003(3) —0.003(3) 

QED - 5.052 -5.088 0.036 

Screened SE 0.381 0.381 

Screened VP, this work - 0.188(2) - 0.188(2) 

Screened VP, Refs. [23, 24] - 0.187(6) - 0.187(6) 

Total -61.320(4) (5) 



Acknowledgments 

Valuable conversations with A. N. Artemyev are gratefully acknowledged. The work re- 
ported in this paper was supported by the Deutsche Forschungsgemeinschaft (Grants No. 
VO 1707/1-1 and PL 254/7-1), by GSI, by RFBR (Grant No. 10-02-00450), by the Russian 
Ministry of Education and Science (Program "Scientific and pedagogical specialists for inno- 
vative Russia", Grant No. P1334), by the grant of the President of the Russian Federation 
(Grant No. MK-3215.2011.2), and by the German- Russian Interdisciplinary Science Center 
(G-RISC) funded by the German Federal Foreign Office via the German Academic Exchange 
Service (DAAD). O.V.A. acknowledges the financial support from the "Dynasty" foundation 
and from the DAAD. The work of D.A.G. was supported by the "Dynasty" foundation and 
by the FAIR - Russia Research Center. O.V.A. and D.A.G. would like to acknowledge the 
hospitality shown by TUD. 



[1] I. Klaft, S. Borneis, T. Engel, B. Fricke, R. Grieser, G. Huber, T. Kiihl, D. Marx, R. Neumann, 
S. Schroder, P. Seelig, and L. Volker, Phys. Rev. Lett. 73, 2425 (1994). 



15 



[2] J. R. Crespo Lopez-Urrutia, P. Beiersdorfer, D. W. Savin, and K. Widmann, Phys. Rev. Lett. 
77, 826 (1996). 

[3] J. R. Crespo Lopez-Urrutia, P. Beiersdorfer, K. Widmann, B. B. Birkett, A.-M. Martensson- 

Pendrill, and M. G. H. Gustavsson, Phys. Rev. A 57, 879 (1998). 
[4] P. Seelig, S. Borneis, A. Dax, T. Engel, S. Faber, M. Gerlach, C. Holbrow, G. Huber, T. Kiihl, 

D. Marx, K. Meier, P. Merz, W. Quint, F. Schmitt, M. Tomaselli, L. Volker, H. Winter, 

M. Wiirtz, K. Beckert, B. Franzke, F. Nolden, H. Reich, M. Steck, and T. Winkler, Phys. 

Rev. Lett. 81, 4824 (1998). 
[5] P. Beiersdorfer, S. B. Utter, K. L. Wong, J. R. Crespo Lopcz-Urrutia, J. A. Britten, H. Chen, 

C. L. Harris, R. S. Thoe, D. B. Thorn, E. Trabert, M. G. H. Gustavsson, C. Forssen, and 

A.-M. Martensson-Pendrill, Phys. Rev. A 64, 032506 (2001). 
[6] S. M. Schneider, W. Greiner, and G. Soff, Phys. Rev. A 50, 118 (1994). 
[7] V. M. Shabaev, J. Phys. B 27, 5825 (1994). 

[8] V. M. Shabaev, M. B. Shabacva, and I. I. Tupitsyn, Phys. Rev. A 52, 3686 (1995). 
[9] H. Persson, S. M. Schneider, W. Greiner, G. Soff, and I. Lindgren, Phys. Rev. Lett. 76, 1433 
(1996). 

[10] S. A. Blundell, K. T. Cheng, and J. Sapirstein, Phys. Rev. A 55, 1857 (1997). 

[11] V. M. Shabaev, M. Tbmaselh, T. Kiihl, A. N. Artemyev, and V. A. Yerokhin, Phys. Rev. A 
56, 252 (1997). 

[12] P. Sunnergren, H. Persson, S. Salomonson, S. M. Schneider, I. Lindgren, and G. Soff, Phys. 
Rev. A 58, 1055 (1998). 

[13] V. M. Shabaev, M. B. Shabaeva, I. I. Tupitsyn, and V. A. Yerokhin, Hyperfine Interact. 14, 
129 (1998). 

[14] S. Boucard and P. Indelicato, Eur. Phys. J. D 8, 59 (2000). 

[15] V. M. Shabaev, A. N. Artemyev, O. M. Zherebtsov, V. A. Yerokhin, G. Plunien, and G. Soff, 

Hyperfine Interact. 27, 279 (2000). 
[16] J. Sapirstein and K. T. Cheng, Phys. Rev. A 63, 032506 (2001). 

[17] V. M. Shabaev, A. N. Artemyev, V. A. Yerokhin, O. M. Zherebtsov, and G. Soff, Phys. Rev. 
Lett. 86, 3959 (2001). 

[18] A. N. Artemyev, V. M. Shabaev, G. Plunien, G. Soff, and V. A. Yerokhin, Phys. Rev. A 63, 
062504 (2001). 

16 



[19] V. A. Yerokhin and V. M. Shabaev, Phys. Rev. A 64, 012506 (2001). 

[20] Z. Andjelkovic, S. Bharadia, B. Sommer, M. Vogel and W. Nortershauser, Hyperfine Interact. 

196, 81 (2010). 
[21] W. Nortershauser, Hyperfine Interact. 199, 131 (2011). 
[22] W. Nortershauser, private communication. 

[23] A. V. Volotka, D. A. Glazov, V. M. Shabaev, I. I. Tupitsyn, and G. Plunien, Phys. Rev. Lett. 
103, 033005 (2009). 

[24] D. A. Glazov, A. V. Volotka, V. M. Shabaev, I. I. Tupitsyn, and G. Plunien, Phys. Rev. A 
81, 062112 (2010). 

[25] A. V. Volotka, D. A. Glazov, O. V. Andreev, V. M. Shabaev, I. I. Tupitsyn, and G. Plunien, 

to be published in Phys. Rev. Lett. 
[26] J. Sapirstein and K. T. Cheng, Phys. Rev. A 67, 022512 (2003); 74, 042513 (2006); 78, 022515 

(2008). 

[27] N. S. Oreshkina, A. V. Volotka, D. A. Glazov, I. I. Tupitsyn, V. M. Shabaev, and G. Plunien, 

Opt. Spektrosk. 102, 889 (2007) [Opt. Spectrosc. 102, 815 (2007)]. 
[28] Y. S. Kozhedub, D. A. Glazov, A. N. Artemyev, N. S. Oreshkina, V. M. Shabaev, 1. 1. Tupitsyn, 

A. V. Volotka, and G. Plunien, Phys. Rev. A 76, 012511 (2007). 
[29] N. S. Oreshkina, D. A. Glazov, A. V. Volotka, V. M. Shabaev, I. I. Tupitsyn, and G. Plunien, 

Phys. Lett. A 372, 675 (2008). 
[30] A. V. Volotka, D. A. Glazov, I. I. Tupitsyn, N. S. Oreshkina, G. Plunien, and V. M. Shabaev, 

Phys. Rev. A 78, 062507 (2008). 
[31] A. V. Volotka, V. M. Shabaev, G. Plunien, and G. Soff Eur. Phys. J. D 23, 51 (2003). 
[32] V. M. Shabaev, Phys. Rep. 356, 119 (2002). 
[33] P. J. Mohr, G. Plunien, and G. Soff, Phys. Rep. 293, 227 (1998). 
[34] M. Gyulassy, Nucl. Phys. A 244, 497 (1975). 
[35] G. A. Rinker and L. Wilets, Phys. Rev. A 12, 748 (1975). 
[36] G. SofF and P. Mohr, Phys. Rev. A 38, 5066 (1988). 

[37] A. N. Artemyev, T. Beier, G. Plunien, V. M. Shabaev, G. Soff, and V. A. Yerokhin, Phys. 
Rev. A 60, 45 (1999). 

[38] A. I. Milstein and A. S. Yelkhovsky, Phys. Lett. B 233, 11 (1989); Zh. Eksp. Teor. Fiz. 99, 
1068 (1991) [Sov. Phys. JETP 72, 592 (1991)]. 

17 



[39] V. M. Shabaev, I. I. Tupitsyn, V. A. Yerokhin, G. Plunien, and G. Soff, Phys. Rev. Lett. 93, 
130405 (2004). 

[40] J. Sapirstein and W. R. Johnson, J. Phys. B 29, 5213 (1996). 

[41] A. G. Fainshtein, N. L. Manakov, and A. A. Nckipelov, J. Phys. B 23, 559 (1990). 

[42] I. Angeh, At. Data Nucl. Data Tables 87, 185 (2004). 

[43] N. J. Stone, At. Data Nucl. Data Tables 90, 75 (2005). 

[44] A. V. Nefiodov, G. Plunien, and G. SofF, Phys. Lett. B 552, 35 (2003). 



18 



